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We present a measurement of the branching fraction and photon energy spectrum for the decay 
B — > Xs7 using data from the BABAR experiment. The data sample corresponds to an integrated 
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luminosity of 210 fb~^, from which approximately 680 000 BB events are tagged by a fully recon- 
structed hadronic decay of one of the B mesons. In the decay of the second B meson, an isolated 
high-energy photon is identified. We measure B{B — > Xs"/) = (3.66 ± 0.85stat ± 0.60syst) x 10"** for 
photon energies E~, above 1.9 GeV in the B rest frame. From the measured spectrum we calculate 
the first and second moments for different minimum photon energies, which are used to extract the 
heavy-quark parameters mb and /ij. In addition, measurements of the direct CP asymmetry and 
isospin asymmetry are presented. 

PACS numbers: 13.20.He, 13.30.Ce, 12.39.Hg 



INTRODUCTION 

We present measurements of the branching fraction 
and photon energy spectrum of the rare radiative pen- 
guin decay B Xg^ using T(4S') BB events. We 
use a new technique where one of the B mesons (called 
the tag B) decays to hadrons and is fully reconstructed. 
This approach allows for the determination of the charge, 
flavor and momentum of both of the B mesons, and thus 
the photon spectrum can be determined in the rest frame 
of the signal B. The method results in an improved pu- 
rity for the signal sample, allows separate measurements 
for charged and neutral B mesons and enables the mea- 
surement of the direct CP asymmetry Acp. This ap- 
proach is complementary to those used in previous stud- 
ies P, 0, H, and incurs different systematic uncertain- 
ties. 

In the Standard Model (SM), the decay b ^ sj pro- 
ceeds via a flavor-changing neutral current. The decay 
is sensitive to new physics through non-SM heavy parti- 
cles entering at the loop level [5|. Recent next-to-next- 
to-leading-order calculations predict SM branching frac- 
tions in the range B{B X,7) = (3.0 - 3.5) x 10"^ 
for > 1.6 GeV with uncertainties that vary from 7% 
to 14% @, 0, i- Here is the energy of the signal 
photon in the rest frame of the B meson, and the cutoff 
is chosen to avoid non-perturbative effects at lower en- 
ergies. The current world average measured branching 
fraction is I3(B^ X^-f) = (3.55 ± 0.26) x 10"^ {E-y > 
1.6 GeV) 0, [l3|. The moments of the photon energy 
spectrum are sensitive to the Heavy Quark Expansion 
parameters mt and /u^, related to the mass and momen- 
tum of the b quark within the B meson [llj. Improved 
measurements of these parameters can be used to reduce 
the uncertainty in the CKM matrix elements |ycfc| and 
\Vub\ 0,0. 

The measurements presented here are based on a sam- 
ple of 232 million BB pairs collected on the T(4S') 



resonance by the BaBAR detector [IJ] at the PEP- 
II asymmetric-energy e~^e~ storage ring operating at 
SLAG, corresponding to an integrated luminosity of 
210 fb~^. After reconstruction of the tag B, the remain- 
ing particles in the event are assigned to the second B 
(the signal B) and events containing a high-energy pho- 
ton are selected. The signal process B X^ at this stage 
is taken to mean events from either b sj or b d'y 



decays; the small contribution from b ^ dj is subtracted 
at the end of the analysis. The sample also includes back- 
ground from continuum (non- _Bi?) events and BB events 
in which the tag B is misreconstructed. These are sub- 
tracted by means of a fit to the beam-energy-subsituted 
mass (defined below) of the tag B. 

The remaining background events, where the photon 
candidate is not from the signal process (e.g., a photon 
from a tt" or 7/ decay), are subtracted using a Monte Garlo 
(MG) model based on EvtGen ^ and GEANT4 [11]. The 
MG predictions arc scaled to data in the low Ej region, 
where the signal contribution is very small. This allows a 
reliable measurement for photon energies E^ > 1.9 GeV. 
Finally, to compare with other experiments and predic- 
tions, the measured rate is extrapolated using theoretical 
models to give the rate for E^ > 1.6 GeV. 

This measurement is currently limited by statistics, 
and furthermore, the dominant systematic errors are of 
the type that should decrease with a larger data sam- 
ple. Therefore, the approach followed here is expected 
to provide an increasingly competitive level of precision 
when applied to the larger data sample currently being 
collected by the BaBAR experiment. 



EVENT SELECTION 

Using 1114 exclusive hadronic decay channels [lit , 
which represent about 5% of the total decay width of 
the B'^ and B'^ mesons, we identify events in which one 
of the two B mesons is fully reconstructed. The kine- 
matic consistency of the tag B candidates is checked with 
two variables, the beam-energy-substituted mass toes = 
^s/4 — 0^, and the energy difference AE = Eb — \fsl'2., 
where s is the total energy squared in the center-of-mass 
(cm.) frame, and Eb and pB are the the cm. energy 
and momentum of the tag B candidate. We require 
|A_E| < 60McV, a window of approximately ±3(t. 

Those particles in the event that are not reconstructed 
as part of the tag B are regarded as coming from the sig- 
nal B. Among these particles we require an isolated pho- 
ton candidate with energy E^ > l.i GeV in the B frame. 
To ensure a well reconstructed photon, we require the 
electromagnetic shower to lie well within the calorime- 
ter acceptance and to satisfy isolation and shower shape 
requirements. 
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The background events consist of non-signal B decays 
and continuum background from uu, dd, ss and cc events. 
The continuum events are suppressed by using a Fisher 
discriminant that combines 12 variables related to the 
different event decay topologies of BB and continuum 
events. These include event-shape variables such as the 
thrust, as well as information on the energy flow relative 
to the direction of the candidate signal photon. 

To discriminate against photons from tt*^ and rj decays, 
we combine the signal candidate photon with any other 
photon in the event associated with the signal B. The 
event is vetoed if the pair's invariant mass is consistent 
with a tt'^ or -ij. Furthermore, the event is rejected if the 
candidate photon combined with a tt^ is consistent with 
a ^ tt'^tt^ decay assuming that the second photon 
from the tt^ decay is lost. 



FIT OF SIGNAL RATES 

The distribution of jties for the selected events has a 
peak around the mass of the B meson, corresponding to 
correctly reconstructed BB events, and a broad back- 
ground component that stems from non-BB and mis- 
reconstructed BB events. The peak is modeled with a 
Crystal Ball (CB) function This contains two pa- 

rameters that correspond to the mean and width of the 
Gaussian core and two additional parameters that de- 
scribe a power-law tail extended to masses below the core 
region. The non-peak background term is described with 



an ARGUS function [17 1 



Applying the selection criteria outlined above yields 
approximately 7 700 events. We divide the event sample 
into 14 intervals of photon energy, each 100 MeV wide, 
spanning the range 1.3 to 2.7 GeV. In each interval, we 
extract the number of peak events with a binned maxi- 
mum likelihood fit to the rriEs distribution. 

The limited size of the data sample means that it is 
not possible to fit all of the parameters related to the 
shape of the CB and ARGUS functions individually in 
separate intervals of photon energy. One expects, how- 
ever, a smooth variation of the shapes as a function of 
E^. To impose this smoothness, a simultaneous fit of the 
mES distributions for all of the photon-energy intervals is 
carried out. The variation of the shape parameters with 
photon energy is described by polynomials, whose orders 
are the lowest possible that allow an adequate modeling 
of the data. Examples of the m-Es distributions and re- 
sults of the simultaneous fit are shown in Fig. [T] The 
global is 330 for the charged B sample and 357 for 
the neutral sample, both for 387 degrees of freedom. 

The measured numbers of B events are shown in 
Fig. [1] (c) as a function of photon energy. The points are 
from data; the solid histogram is from a BB MC sample 
that excludes the signal decay B — ^ Xj. Due to the large 
background at low energy the signal region is defined as 



Ej > 1.9 GcV. This choice was optimized in MC studies. 
The MC prediction has been scaled by fitting to the data 
region between 1.3 < E^ < 1.9 GeV, taking into account 
the small contribution from B —f X"f decays in that re- 
gion. For E^ > 1.9 GeV, we observe 119 ±22 B ^ X-f 
signal events over a BB background of 145 ± 9 events. 

For 1.3 < Ej < 1.9 GeV a comparison of the data and 
background gives a of 9.7 for five degrees of freedom. 
The probability to observe a value at least this great is 
8.4%. Our estimate of the systematic uncertainty in the 
background (described below) is in fact smaller than the 
observed data-background difference; therefore we regard 
this difference primarily as a statistical fluctuation. 

To determine the partial branching fractions, we re- 
quire the total number of BB events in the sample after 
selection of the tag B candidates. In a procedure anal- 
ogous to that described for the toes fits in bins of E^, 
we divide the data into four intervals of estimated tag 
B candidate purity and perform a simultaneous fit of the 
Toes distributions. We obtain approximately 680 000 BB 
events corresponding to an efficiency of 0.3%. 



DETERMINING THE PHOTON SPECTRUM 



The differential decay rate {l/rB){dT/dEj) is mea- 
sured in bins of the (B-frame) photon energy for Ej > 
1.9 GeV up to the kinematic limit at 2.6 GeV. It is esti- 
mated for the ith bin as 



1 dV, 

Tb dE^ 



s^Nb 



(1) 



where Ni is the number of B events in the bin, hi is the 
number of B mesons from decays other than B — > X7, 
Nb is the total number of B mesons in the sample, and 
Ei is the efficiency, which corrects for both acceptance 
and bin-to-bin resolution effects. The values bi are de- 
termined by means of a simultaneous fit to the TOes dis- 
tributions as described previously, using a sample of MC 
data consisting of BB events excluding the signal decay 
B Xj. As the differential decay rate is normalized 
using the total width of the B meson, Tb, the integral 
of ^ over all photon energies yields the branching frac- 
tion. To evaluate the selection efficiency e^, we model the 
signal photon energy spectrum using the kinetic scheme 



[18j with TOh = 4.60 GcV and 
of S j is determined from 



0.4 GcV^ The value 



N. 



found,'! /-^sim 



tag , 



(2) 



gon 



where iVfound,i is the number of events found in a MC 
sample oi B Xsj with detector simulation and Ngim 
is the number of events in the simulated sample. These 
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FIG. 1: Fits to the distribution of the beam-energy-substituted mass mEs for two E.y regions. The dashed curve shows the 
CB term and the dotted curve is the ARGUS term, corresponding to B and non-B events, respectively; the sohd curve is their 
sum. (a) 1.6 GeV < Ej < 1.7 GeV for the charged B sample, (b) 2.3 GeV < < 2.4 GeV for the neutral B sample. (c_)_ The 
measured numbers of B events as a function of photon energy. The points are from data; the histogram is from a BB MC 
sample which excludes the signal decay B — > X^. 



quantities are found using the same fit procedure as ap- 
plied to the real data for Ni and Nb- In the denomi- 
nator of -/VtruG,i is the true number of events with 
photon energies in bin i and iVgen is the total number 
of events generated. These values are determined us- 
ing the event generator for B Xs'^ decays only, with- 
out detector simulation. The factor Ctag, estimated us- 
ing the MC model, corrects for the small dependence 
of the probability to find a tag B on the presence of a 
B X7 final state. The efficiency increases roughly 
linearly with photon energy, and is approximately 30% 
(65%) for = 1.9 GeV (2.6 GeV). 

To compare with other results we subtract the B — > 
Xd"f component from the differential decay rates using 
the Standard Model prediction (for the CP and isospin 
asymmetries discussed below, however, we do not make 
this correction). The values B{B Xd^) and B{B 
Xsj) are in the ratio |t4d/VtsP assuming the same effi- 
ciency for the two categories of events. Therefore, the 
branching ratio is lowered by (4.0 ± 0.4)% [11, [lO]. 

SYSTEMATIC UNCERTAINTIES 

There are four main sources of systematic uncertainty, 
which are summarized in Table [H modeling of the BB 
background, the uies fits, detector response and depen- 
dence on the B ^ Xs'j signal model. In addition there 
is an uncertainty from the subtraction of the B — > X^j 
contribution. 

After subtraction of the non-peak background us- 
ing the TTiEs distribution, the remaining background is 
mainly composed of BB events with the selected photon 
coming from a 7r° or 77 decay. Photons from tt" account 
for 55% to 65% depending on and the charge of the tag 
B, while the contribution from rj mesons varies from 18% 
to 29%. The remaining backgrounds include fake photons 
from h annihilation, real photons from bremsstrahlung 



or from u decays, and electromagnetic showers from e 
misidentified as photons. As the MC prediction for the 
BB background is scaled to the data at low energy, there 
is no uncertainty stemming from the absolute rate, but 
rather only from the shape of the distribution as a func- 
tion of . The uncertainty from the inclusive tt" and 77 
spectra is investigated by using E-y dependent correction 
factors for the 7r° and 77 yields from a large control sam- 
ple oi B Xj candidate events, obtained using a lepton 
tag. These correction factors are typically around 5% for 
7r° yields while they can be up to 30% for 77 yields. The 
remaining backgrounds have a roughly linear slope with 
E^; this is varied by ±30%. We use the difference ob- 
tained with the modified MC compared to the standard 
MC simulation as a systematic uncertainty. 

To assess the uncertainty related to the parameteri- 
zation chosen for the tties fit, additional coefficients are 
introduced that allow linear or higher-order dependence 
of the CB and ARGUS function shape parameters on 
the photon energy. The maximum variation in the fitted 
rates is taken as the systematic uncertainty. A similar set 
of variations for the dependence of the shape parameters 
on the B meson purity is carried out for the 77t,es fits used 
to determine the total number of B mesons in the data 
sample. To allow for a small peaking component in the 
distribution of 77t,es from B^ decays reconstructed as B^ 
(-B") decays and vice versa, we remove these events from 
the MC sample and take the difference in the result as a 
systematic uncertainty. 

The uncertainties related to the detector modeling and 
event reconstruction are estimated by comparing MC 
simulations of track and photon efficiencies as well as 
particle identification efficiencies with data control sam- 
ples. From these comparisons we estimate corresponding 
systematic errors, which are in all cases small compared 
to other uncertainties. 

To assess the uncertainty in the efficiency due to the 
assumed shape of the E^ spectrum, we vary mb and 
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FIG. 2: The partial branching fractions {l/rB)idr/dEj) 
with statistical (inner) and total (outer) uncertainties. 



in the kinetic scheme by ±0.1 GcV and ±0.1 GeV , re- 
spectively. These variations are large compared to the 
uncertainties in the world average [lO| in order to cover 
alternative Ansatze for the the heavy quark distribution 
function [2l|, . They also account for uncertainties re- 
lated to the small rate of i? ^ X'y decays expected below 
1.9 GcV. 



TABLE I: Results for the differential decay rate 
{l/TB){dr/dE^) and moments of the photon spectrum 
with statistical and systematic errors. The major contri- 
butions to the systematic uncertainties are also listed: (a) 
background modeling, (b) tties fit parameterization, (c) 
detector response, (d) B — > Xs^ model. 

(i/rfl)(dr/d£;^)(io-") 



E-y ( GeV) Value 


Cstat 


•^syst 


(a) 


(b) 


(c) 


(d) 


1.9-2.0 


0.28 


0.56 


0.34 


0.26 


0.13 


0.19 


0.03 


2.0-2.1 


0.60 


0.42 


0.24 


0.18 


0.12 


0.08 


0.05 


2.1-2.2 


0.31 


0.29 


0.14 


0.11 


0.06 


0.03 


0.03 


2.2-2.3 


0.40 


0.23 


0.13 


0.07 


0.05 


0.09 


0.03 


2.3-2.4 


0.91 


0.22 


0.13 


0.07 


0.08 


0.05 


0.06 


2.4-2.5 


0.74 


0.17 


0.09 


0.05 


0.05 


0.02 


0.05 


2.5-2.6 


0.43 


0.12 


0.09 


0.03 


0.03 


0.07 


0.04 








(E^) (GeV) 








E-y ( GeV) Value 


Cstat 


f^syst 


(a) 


(b) 


(c) 


(d) 


1.9-2.6 


2.289 


0.058 


0.027 


0.018 


0.019 


0.009 


0.002 


2.0-2.6 


2.315 


0.036 


0.019 


0.013 


0.011 


0.009 


0.001 


2.1-2.6 


2.371 


0.025 


0.009 


0.007 


0.005 


0.003 


0.001 


2.2-2.6 


2.398 


0.016 


0.004 


0.003 


0.003 


0.001 


0.000 


2.3-2.6 


2.427 


0.010 


0.006 


0.000 


0.001 


0.005 


0.000 








-{E,)r) (GeV^) 






E-, ( GeV) Value 


Cstat 


*7syst 


(a) 


(b) 


(c) 


(d) 


1.9-2.6 


0.0334 


0.0124 


0.0062 


0.0040 


0.0025 


0.0037 


0.0013 


2.0-2.6 


0.0265 


0.0057 


0.0024 0.0018 


0.0010 


0.0007 


0.0011 


2.1-2.6 


0.0142 


0.0037 


0.0013 


0.0009 


0.0005 


0.0004 


0.0006 


2.2-2.6 


0.0092 


0.0015 


0.0010 


0.0002 


0.0002 


0.0009 


0.0003 


2.3-2.6 


0.0059 


0.0007 


0.0003 


0.0000 


0.0000 


0.0003 


0.0002 



RESULTS 



The partial branching fractions {l/rB){dT/dE^) are 
shown in Fig. [2] after all corrections. The inner error bars 
show the statistical uncertainties. The outer error bars 
show the quadratic sum of the statistical and systematic 
terms. By integrating the spectrum, we obtain B{B 
Xs-/) = (3.66 ± 0.85stat ± O.eOsyst) x lO""^. The results 
for the differential decay rate and for the moments of 
the photon energy spectrum for various minimum photon 
energies £'cut are given in TableHl The branching fraction 
for larger values of -Ecut and the correlations between the 
measurements are given in Tables IlllfVl Our results are in 
good agreement with those presented in Refs. [H, 

We also measure the isospin asymmetry Aq-, 



Ao- = 



r(BO ^ XsA 7) - r(s- ^ XsA 7) 



r(B" ^ x,^d 7) + r(s- ^ x,^d 7) 



(3) 



where inclusion of charge conjugate modes is implied. It 
has been argued that enhanced power corrections to the 
B — > Xsj rate could also lead to values of Aq- as large 
as +10% 23|. Therefore, experimental measurements of 
Aq- can help determine the size of these effects and hence 
reduce the theoretical uncertainty on the total rate. To 
obtain decay rates from the branching fractions we use 
the B meson hfetimes: t{B°) = 1.530 ± 0.008 ps and 
t{B+ ) = 1.638 ± 0.011 ps For photon energies 



greater than 2.2 GeV, we obtain Aq- = — 0.06±0.15stat± 

0.07syst. 

The direct CP asymmetry Acp, 

^ ^ B{B^X,,dj)-B{B^XsMl) 1 
'''' 6(B^X,,d7)+S(S^X3,d7)l-2w' 

is measured by splitting the tag sample into B and B 
mesons. The dilution factor 



- accounts for the mistag 
fraction w, here simply the time integrated B'-' mixing 
probability of Xd = 0.188 ± 0.003 [H multiplied by the 
fraction of B^ events in the total data sample. Aqp can 
be significantly enhanced by new physics Il9| while in the 
SM it is predicted to be around 10~^ [2^,l26(. We obtain 
a value of Acp = 0.10 ± O.lSstat ± 0.05syst for photon 
energies above 2.2 GeV. 

For both Aq- and Acp, a photon energy cutoff of 
2.2 GeV is chosen because it facilitates comparison with 
previous results and minimizes the total uncertainty. Our 
results are in goo d ag reement with previous measure- 



Mm 



Finally, we use heavy quark expansions in the kinetic 
scheme [l8| and our measurements of the moments 
to determine the parameters mi, and We include 
the theoretical uncertainties quoted in Ref. [l3| in the 
overall covariancc matrix used in the fit. To minimize 
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the theoretical uncertainty we only use moments with 
Ecut < 2.0 GeV and obtain rub = 4.46+1^:23 GeV and /i^ = 
0.64+°J^ GeV^ with a correlation of p = -0.94. 



CONCLUSIONS 

We have measured the B X^j branching fraction 
and moments of the photon energy spectrum above sev- 
eral minimum photon energies. We find B{B Xsj) = 
(3.66 ± 0.85stat ± O.eOsyst) X 10"^ for photon energies 
above 1.9 GeV. Dividing by an extrapolation fac- 
tor of 0.936 ± 0.010 we obtain B{B Xs-f) = 
(3.91±0.91stat±0.64syst) X 10-4 for > 1.6 GeV. The 
moments of the spectrum can be used to improve the 
knowledge of the heavy quark parameters mf, and ; we 

+0.21, 



obtain mf, = 4.46z:o;23 GeV and 



0.64t°;i^^ GeV^ in 
the kinetic scheme. In addition we measured the isospin 
asymmetry Aq- = —0.06 ± O.lSstat i O.OTgyst and di- 
rect CP asymmetry Acp = 0.10 ± O.lSstat ± 0.05syst for 
photon energies above 2.2 GcV. The full reconstruction 
(recoil) method provides an almost background free mea- 
surement above photon energies of 2.2 GeV. Although 
statistics are limited at present, this approach is ex- 
pected to provide a competitive measurement of the de- 
cay B — > with the larger data sample that is be- 
ing accumulated at the B -Factories, in particular as the 
main systematic uncertainties will also be reduced with 
a larger data sample. 
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TABLE 11: Results for B{B Xs'y) for different minimum photon energies with statistical and systematic errors. Details on 
the major contributions to the systematic uncertainties are also given. 



B{B ^ Xs'y) (10-*) 



-E-y range 


Value 


''"stat 




Background 


mES fit 


Detector B 




(GeV) 








modeling 


parameterization response 


model 


1.9-2.6 


3.66 


0.85 


0.60 


0.35 


0.45 


0.18 


0.08 


2.0-2.6 


3.39 


0.64 


0.47 


0.31 


0.34 


0.07 


0.06 


2.1-2.6 


2.78 


0.48 


0.35 


0.22 


0.24 


0.08 


0.05 


2.2-2.6 


2.48 


0.38 


0.27 


0.14 


0.19 


0.10 


0.05 


2.3-2.6 


2.07 


0.30 


0.20 


0.10 


0.15 


0.04 


0.05 



TABLE IIL Correlations between the systematic uncertainties for the differential decay rate measurements. (The statistical 
correlations between Ery bins are negligible.) 



E.y interval 






(1/r 


B){dr/dE^) 






(GeV) 


1.9-2.0 


2.0-2.1 


2.1-2.2 


2.2-2.3 2.3-2.4 


2.4-2.5 


2.5-2.6 


1.9-2.0 


1.000 


-0.004 


0.311 


0.949 -0.108 


0.083 


0.849 


2.0-2.1 




1.000 


0.912 


0.096 0.805 


0.721 


-0.087 


2.1-2.2 






1.000 


0.352 0.712 


0.699 


0.152 


2.2-2.3 








1.000 0.111 


0.310 


0.940 


2.3-2.4 








1.000 


0.969 


0.115 


2.4-2.5 










1.000 


0.341 


2.5-2.6 












1.000 



TABLE IV: Statistical correlations between (Ej) and {{Ej — (Ej))'^) measurements with different minimum photon energies E. 



Ej range 
(GoV) 


1.9-2.6 2.0-2.6 2.1-2.6 2.2-2.6 2.3-2.6 


{{E-, - {E,)Y) 
1.9-2.6 2.0-2.6 2.1-2.6 2.2-2.6 2.3-2.6 


1.9-2.6 
2.0-2.6 
{E^} 2.1-2.6 

2.2- 2.6 

2.3- 2.6 


1.000 0.503 0.195 0.053 -0.001 
1.000 0.418 0.141 0.027 
1.000 0.427 0.157 
1.000 0.408 
1.000 


-0.897 -0.419 -0.180 -0.040 0.041 
-0.310 -0.807 -0.368 -0.093 0.069 
-0.124 -0.342 -0.822 -0.244 0.109 
-0.017 -0.054 -0.153 -0.550 0.200 
0.008 0.019 0.032 0.095 0.363 


1.9-2.6 
2.0-2.6 

{{E^ - {E^)f) 2.1-2.6 

2.2- 2.6 

2.3- 2.6 




1.000 0.266 -0.041 -0.048 0.054 
1.000 0.094 -0.052 0.113 
1.000 0.177 0.144 
1.000 0.350 
1.000 



TABLE V: Systematic correlations between {E-y) and {{E-y — {Ej))^) measurements with different minimum photon energies E. 



E~f range 
(GeV) 


1.9-2.6 2.0-2.6 2.1-2.6 2.2-2.6 2.3-2.6 


{{E-, - {E,)Y) 
1.9-2.6 2.0-2.6 2.1-2.6 2.2-2.6 2.3-2.6 


1.9-2.6 
2.0-2.6 
(E^) 2.1-2.6 

2.2- 2.6 

2.3- 2.6 


1.000 0.556 0.811 0.707 0.128 
1.000 0.946 0.667 0.751 
1.000 0.765 0.653 
1.000 0.527 
1.000 


0.298 0.875 0.599 -0.030 -0.013 
0.903 0.511 0.943 0.775 0.724 
0.741 0.771 0.932 0.553 0.546 
0.412 0.777 0.731 0.307 0.416 
0.745 0.102 0.658 0.911 0.964 


1.9-2.6 
2.0-2.6 

{(E^ - {E,,)f) 2.1-2.6 

2.2- 2.6 

2.3- 2.6 




1.000 0.253 0.841 0.868 0.770 
1.000 0.656 -0.059 -0.019 
1.000 0.627 0.597 
1.000 0.967 
1.000 



